Anderson impurity model at finite Coulomb interaction U: generalized Non-crossing 

Approximation 
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We present an extension of the non-crossing approximation (NCA), which is widely used to 
calculate properties of Anderson impurity models in the limit of infinite Coulomb repulsion U — > oo, 
to the case of finite U. A self-consistent conserving pseudo-particle representation is derived by 
symmetrizing the usual NCA diagrams with respect to empty and doubly occupied local states. 
This requires an infinite summation of skeleton diagrams in the generating functional thus defining 
the "Symmetrized finite-U NCA" (SUNCA). We show that within SUNCA the low energy scale T K 
(Kondo temperature) is correctly obtained, in contrast to other simpler approximations discussed 
in the literature. 
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I. INTRODUCTION 

Anderson impurity models have been of considerable 
interest recently as generic models of local systems with 
internal degrees of freedom coupled to a Fermi gas. Al- 
though fiest introduced as models for magnetic impurities 
in metalad, they describe two-level systems in metalsH, 
quantum dots in mesoscopic structuresoQ and strongly 
correlated lattice systems in the local approximation of 
the Dynamic Mean Field Theory (DMFT)eI as well. In a 
nutshell, the Anderson model features one or several lo- 
cal levels hybridizing with the conduction electron states 
of the metal. Multiple occupancy of the local levels is 
inhibited by the strong Coulomb repulsion U between 
electrons in the local states. As a consequence, the local 
levels are approximately singly occupied, giving rise to 
a magnetic moment or an equivalent degree of freedom. 
Due to the coupling to the conduction electron system, 
the local moment is screenedEl, or in a multi-channel sit- 
uation forms a more complicated many-body statecl. 

Most investigations of Anderson models have concen- 
trated on the case of infinite Coulomb repulsion U. The 
corresponding restriction of the local Hilbert space to 
electron occupation n e i < 1 allows for an economical 
treatment in terms of pseudo-particle representationsQIa 
and a projection onto the physical sector of Hilbert space. 
In this framework, the simplest approximation consists 
of second-order self-consistent perturbation theory in the 
hybridization, the so-called Non-Crossing Approximation 
(NCApil Although the NCA has its limitations, it is 
a valuable tool for extracting the complex many-body 
physics of Anderson impurity models. In the single chan- 
nel case the NCA accounts correctly for the formation of 
a Kondo resonance at the Fermi level below the Kondo 
temperature Tr-EJ, even though the appearance of a local 
Fermi liquid state at temperatures T <C Tk is not cap- 
tured in this approximation! 2 ]. In the multichannel case 
even the correct low temperature power law behavior is 
obtained in NCAE3. However, in order to capture, e.g., 



the physics of the upper and the lower Hubbard bands 
in a DMFT description of the Hubbard model and the 
Mott-Hubbard metal-insulator transition, it is essential 
to consider the case of large but finite U. It is therefore 
desirable to develop a generalization of NCA to the case 
of finite Coulomb interaction. In the following we present 
a straightforward generalization of NCA, which conserves 
the symmetry of virtual transitions to the empty local 
level or doubly occupied local level states. This is essen- 
tial for recovering the correct Kondo [temperature Tk, 
as pointed out by Pruschke and Grewai-2 and, as will be 
shown, requires an infinite summation of a certain class of 
crossing diagrams. We find that inclusion of only the first 
crossing term in this resummatioro, while contributing 
the larger part of the change of Tk, is not sufficient to 
provide a qualitatively correct Kondo temperature. 



II. PSEUDOPARTICLE REPRESENTATION OF 
THE MODEL 



The model we consider describes a local impurity level 
(called d-level in the following), hybridizing with the con- 
duction electron states. The energy Ed of the level may 
be located below or above the Fermi energy. Two elec- 
trons with spins f and J, on the local level experience a 
Coulomb interaction U. The local states will be assumed 
to be created by pseudo-particle operators fl (singly oc- 
cupied state with spin a), b* (empty state) and (doubly 
occupied state) acting on a vacuum state without any im- 
purity. We choose f a to be fermion and a, b to be boson 
operators, where b will be called the "light" and a the 
"heavy" boson. The creation operator for the local phys- 
ical electron can then be written as S a — /J6 + era* f— a , 
where the pseudo-particle occupation numbers must sat- 
isfy the operator constraint 
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expressing the fact that at any instant of time the im- 
purity is in exactly one charge state, empty, singly, or 
doubly occupied, respectively. The fermion operators ct 

create electrons in conduction electron states | ka > with 
energy e^. The Hamiltonian then takes the form 



+ Ua^a + V ( c \/f° + ac lj-* a + h - c ) ' ( 2 ) 



where V is the hybridization matrix element. For later 
use we define the conduction electron density of states at 
the Fermi energy as 7V(0) and the effective hybridization 
T = ttN{Q)V 2 . 



III. GAUGE SYMMETRY AND PROJECTION 
ONTO THE PHYSICAL HILBERT SPACE 

The model described by the auxiliary particle Hamil- 
tonian (g) is invariant under simultaneous, local U{\) 
gauge transformations, f a — > /o-e^W, b — > be 1 ^ 1 ', a — ► 
ae l ^ l \ where 4>{t) is an arbitrary, time-dependent phase. 
This gauge symmetry guarantees the conservation of the 
local charge Q in time. In order to project onto the phys- 
ical subspace Q — 1 , it is therefore sufficient to carry out 
the projection at time t — > — oo, if the gauge symmetry is 
implemented exactly. One starts with the grandcanonical 
ensemble with respect to Q and the associated chemical 
potential —A. The projection is achieved by taking the 
limit A — > oo of any grandcanonical expectation value 
of a physical operator A acting in the impurity Hilbert 
space 



(A) 
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Here the subscript G denotes the grandcanonical ensem- 
ble. The extra factor Q in the denominator of Eq. (|J) has 
been introduced to project out the Q — subspace. Note 
that in the numerator this factor can be omitted, since 
any physical operator A acting on the impurity states 
consists of powers of d\, d a , which annihilate any state 
in the Q = subspace, d\\Q = 0) = 0, d a \Q = 0) = 0. A 
detailepLdescription of the projection procedure is given 
in RefJiJ. Expectation values in the grandcanonical en- 
semble may be calculated straightforwardly in perturba- 
tion theory in the hybridization V, making use of Wick's 
theorem. The usual resummation techniques may be ap- 
plied. Thus the imaginary time single particle Green's 
functions 



G /CT (T 1 -r 2 ) = -(f[/ CT (r 1 )/t(r 2 )]) ( 



(4) 



and analagously for the two bosons a, b, may be expressed 
in terms of the self-energies Yif^ c {iu>) as 



Gf a (iu)) = [w - A - E d - £/(«*>)] 1 
Gf,(i(jj) = [iw — A — Et,(i^)] _1 (5) 
G„(iw) = [iw - A - 2E d - U - EaMp 1 . 

The local conduction electron Green's function is given 

by 



Gr 



with 



(6) 



(7) 



The physical d-electron Green's function is propor- 
tional to the single-particle conduction electron i-matrix 
t ca (iuj), and is related to the grandcanonical (unpro- 
jected) E ccrG as 

G da (iuj) = -^t ca (iuj) = i \im e /5A E C(T G (iuj, A) , (8) 



V 2 A 



where (3 is the inverse temperature. The physical (pro- 
jected onto the Q — 1 subspace) local conduction electron 
self-energy is then obtained from the t-matrix as 



V 2 G da {iuj) 



l + V 2 G* a {iuj)G da {iLu) ' 



(9) 



IV. GENERATING FUNCTIONAL 

Gauge invariant approximations conserving the local 
charge Q may be derived from a Luttinger-Ward gener- 
ating functional For a given approximation the func- 
tional $ is defined by a sum of closed skeleton diagrams. 
The self-energies E M , /i = a, b, /, c, are obtained by tak- 
ing the functional derivatives 



6$ 
SG~. 



(10) 



The "Non-crossing Approximation" (NCA) in the limit 
U — > oo is defined by the single lowest order diagram 
(2nd order in V) containing a light boson line (the first 
diagram of Fig. 1). In the limit of small hybridization 
element Vr, it appears to be justified to keep only the 
lowest QxcLei: contribution in V. However, as discussed 
in RefeMZl, the singular behavior of vertex functions 
may require to include these as well. This turns out to 
be necessary in the single channel model where the for- 
mation of a many-body resonance state is essential for 
recovering the Fermi liquid behavior, and less so in the 
multi-channel models. Including an infinite class of skele- 
ton diagrams in $ (in a "Conserving T-matrix Approx- 
imation" (CTMA)), which allows to capture a singular 
structure in the spin and charge excitation sectors, the 
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low temperature Fermi liquid-phase of the single channel 
Anderson model is recovered^. 

Here we are interested in constructing a simpler gen- 
eralization of NCA to describe the case of finite U. It 
seems straightforward to define such an approximation 
on the NCA level by adding to the second order skeleton 
diagram for $ containing the light boson (the first dia- 
gram in Fig. la)) the corresponding diagram containing 
the heavy boson (the second diagram in Fig. 1 a)). This 
approximation and certain extensions motivated by per- 
turbative argumentsEiLpr by a 1 /N expansion (N beeing 
the spin degeneracy)£3 have been considered sometime 
ago. However, in the case of finite U NCA was found 
to fail badly: Not even the Kondo energy scale is recov- 
ered in the so-defined approximation. The reason for this 

failure is obvious: In the Kondo regime (n<j<l) the local 
spin is coupled to the conduction electron spin density 
at the impurity through the antiferromagnetic exchange 
coupling 

The two terms on the r.h.s. of this relation arise from 
virtual transitions into the empty and doubly occupied 
local level, which e.g. contribute equally in the symmet- 
ric case \ Ed \= Ed + U. The symmetric occurrence of 
these two virtual processes in all intermediate states is 
not included in the simple extension of NCA proposed 
above. Rather, the self-energy insertions in each of the 
two diagrams contain always only one of the processes, 
leading to an effective J which is only one half of the 
correct value. Correspondingly, the Kondo temperature 
T K ~ exp[— l/(2A/"(0) J)] comes out to scale as the square 
of the correct value, which can be orders of magnitude 
too small. 

To correct this deficiency it is necessary to include ad- 
ditional diagrams, restoring the symmetry between the 
two virtual processes. As a first step one may add the 
next order skeleton diagram to $ (Fig. 0b)). As we 
will show below, this approximation, later referred to as 
UNCA, helps to recover a large part of the correct behav- 
ior of TrO. However, as seen from the preceding discus- 
sion, for a completely symmetric treatment of empty and 
doubly occupied intermediate states one must first con- 
sider the diagrams of bare perturbation theory instead 
of skeleton diagrams. A symmetric class of diagrams is 
generated by replacing a light boson line with a heavy 
boson line in each of the bare (non-skeleton) diagrams 
comprising the NCA, and vice versa. Each replacement 
leads to a crossing of conduction electron lines spanning 
one fermion and at most two boson lines. A conserving 
approximation is then constructed by substituting renor- 
malized propagators for the bare ones and keeping only 
skeleton diagrams. The resulting generating functional 
$ is shown diagrammatically in Fig. [l] a)-c). 




FIG. 1. Diagrammatic representation of the generating 
functional to describe the Anderson impurity model at finite 
U . Throughout this paper, solid, dashed, wiggly and zig-zag 
lines correspond to conduction electron c, pseudo-fermion /, 
light boson b, and heavy boson a propagators, respectively, a) 
NCA including light and heavy boson lines. a)-b) Finite-(7 
NCA (UNCA). This approximation amounts to renormalizing 
only one of the vertices in each of the self-energy diagrams of 
Fig. ^ and keeping only one (light or heavy boson) rung in the 
corresponding vertex function (see text). a)-c) Symmetrized 
finite-f/ NCA (SUNCA). 



These diagrams look similar to the CTMA diagrams men- 
tioned above, but contain one light boson line and an 
arbitrary number of heavy boson lines, or vice versa. 
Diagrams with, for example, two light boson lines and 
an arbitrary number of heavy boson lines (and conduc- 
tion electron lines spanning at most one fermion line) 
are reducible and do not appear. We will call the ap- 
proximation defined by the generating functional given 
by the sum of the diagrams of Fig |l| "Symmetrized 
finite-U NCA" (SUNCA). The above approximation cor- 
responding to the CTMA at U — > oo, termed "Sym- 
metrized finite-U Conserving T-matrix Approximation" 
(SUCTMA) is thus defined in a natural way by summing 
up all skeleton <j> diagrams containing a single closed ring 
of auxiliary particle propagators with an arbitrary num- 
ber of light or heavy boson lines, dressed by conduction 
electron lines spanning only one fermion line. Thus, the 
SUCTMA is defined by adding to the diagrams of the 
SUNCA the CTMA diagrams with (only) light boson 
lines or (only) heavy boson lines. The SUCTMA equa- 
tions have not yet been evaluated. 



3 



V. RESULTS OF SUNCA 



As discussed above, the self-energies E^ are obtained 
by functional differentiation of the generating functional 
with respect to the Green's functions G M . The functional 
<& defined by Fig. |l| leads to an infinite series of diagrams 
for E a , which may be conveniently presented in terms of 
three-point vertex functions (the filled semicircles with 
three legs: one boson and two fermion lines), sec Fig. 
H It is necessary to substract a diagram of 4th order in 
V in each case to avoid double counting. On the level 
of SUNCA and SUCTMA the vertex functions consist of 
ladder summations, with light or heavy boson lines as 
rungs, and are defined diagrammatically in Fig. ||[ Note 
that keeping only a single light or heavy boson rung in 
these vertex functions corresponds to UNCA. 



Eh = 



E„ = 



E f = 



V'Gh = 



+ 



FIG. 2. Diagrammatic representation of the auxiliary par- 
ticle self-energies of SUNCA in terms of the renormalized hy- 
bridization vertices, defined in Fig. ^. In each line the third 
diagram is subtracted in order to avoid double counting of 
terms within the first two diagrams. 

The expressions for the self-energies E M defined by 
Figs. H and [| together with the definition of the Green's 
functions, Eqs. (|5|), constitute a set of nonlinear integral 
equations for E A1 (a;),/i = a,b,f. The local conduction 
electron self-energy E c does not appear in any internal 
Green's functions because it contains at least one aux- 
iliary particle loop, i.e. carries a factor exp(— (3X) and 
thus 



a) 



b) 



+ 



+ 



FIG. 3. Diagrammatic representation of 

the Bethe-Salpeter equations for a) the renormalized light 
boson (empty impurity) and b) the renormalized heavy bo- 
son (doubly occupied impurity) vertex, as generated by the 
SUNCA Luttinger-Ward functional (Fig. 1). 



vanishes due to the exact jpojection onto the physi- 
cal Hilbert space (A — > oojlid. Gd and therefore E c 
may be calculated at the end via Eqs. (0), (||) by us- 
ing the self-consistently determined G a ,a = a,b,f. The 
SUNCA equations are given explicitly in appendix A. Al- 
though these equations are more involved than the regu- 
lar NCA, they are numerically considerably more easily 
tractable than the CTMA equationsBO, since SUNCA 
contains only renormalized three-point vertices (see Fig. 
U) as compared to four-point vertex functions occuring in 
CTMAEI. We solved the SUNCA equations numerically 
in the real frequency representation, i.e. after analytic 
continuation from Matsubara frequencies uj„ to the real 
frequency axis. 

As a first important characteristic feature of the 
pseudo-particles we note that the single-particle excita- 
tion spectrum is powerlaw divergent, G M (w) ~ u>~ a » , = 
a,b, f reflecting the abundance of low energy excitations 
forced by the constraint. At finite temperature T these 
singularities are^ciit off at the scale u> ~ T. As observed 
in earlier worl£3 : E£l, the values of the exponents a M are 
characteristic of the state of the system. In the single 
channel case, when the ground state is a local Fermi liq- 
uid, the exponents may be inferred from the Friedel sum 




FIG. 4. Infrared threshold exponents of the auxiliary par- 
ticle spectral functions —\mGf ab (uj)/-K, in dependence of the 
impurity occupation rid , for fixed values of F and U and vary- 
ing Ed- Dashed curved lines: exact results (Eq. (12)); horizon- 
tal lines: NCA results; data points with error bars: SUNCA 
results. In the Kondo limit (rid — > 1) the exact exponents 
are recovered, while in the mixed valence and empty impu- 
rity regime the SUNCA results for «/ and ab cross over to 
the NCA values. 



rule relating the scattering phase shifts r^ IjCT to the num- 
ber of electrons bound to the impurity in each channel 
A? 



'fJ>,<7 



V/j.a-/^- The exponents in turn are related to 
the T)^^ by the geaeral result first derived for the x-ray 
edge singularities!^, a M = 1 — X^fe^/ 71 ") 2 - This is so, 
since e.g. the heavy boson Green's function, describing 
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the transition amplitude for a doubly occupied impurity 
to be created at time t — and removed at a later time 
is proportional to the overlap of the free Fermi sea of 
conduction electrons with the ground state of the Ander- 
son model. The change in the occupation of the local 
level due to the hybridization with the conduction band 
An MCT is the difference between the t = initial impurity 
occupation n^it = 0) (without hybridization) and the 
occupation in the ground state of the Anderson impurity 
model nda — rid/2, i.e. it depends on the initial condi- 
tions of the different Green's functions G^, fi = a, 6, /. 



Thus we have A 



n,. 



_ 2-n d 



An 



(To is the spin of the fermion in the Green's function 
Gf ao ), and Anb !CT = — The infrared threshold ex- 
ponents of G^w) are therefore given by 



5a,, 



hrid (where 



a a = -1 + 2n d - -4 



for the exponent a a trace the exact behavior in reason- 
able agreement. Clearly the SUNCA does much better 
than the simple NCA. From our experienced with the 
Anderson model in the limit U — > oo, we expect that the 
correct exponents should be recovered in SUCTMA. 

We now turn to the d-electron spectral function 
Ad(ijj) = ^ImGd(w — iO). Fig. || shows the results 
for Ad(uj) for the symmetric Anderson model in the 
Kondo regime (n d ~ 1) at a very low temperature of 
T ~ 10 -2 Tjf. Shown are the results obtained from the 
simple NCA (diagrams of first line in Fig. 1), the pertur- 
batively corrected version UNCA (including the diagram 
in the second line of Fig. 1) and the full SUNCA. The 
inset shows that the width of the Kondo resonance peak, 
which is a measure of Tk, comes out orders of magnitude 
different in the three approximations. In order to com- 
pare the numerical results with the exact expression for 
Tk, 



Oib 



(12) 



In Fig. H we show the exponents a M for different n d , as 
obtained from a numerical solution of the SUNCA equa- 
tions. Also shown is the exact result given by Eq. ( |l2|) 
(dashed lines), and the analytical result that can be ex- 
tracted analytically from the NCA equations (defined 
the first two diagrams in Fig. 1) in analogy to RefJl3. 
The numerical results of SUNCA (data points) are seen 
to approach the exact result in the limit rid — > 1, but in 
the case of ctb and a/ appear to follow the NCA result 
rather than the exact result for < 0.8. The results 
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FIG. 5. Local electron spectral function calculated using 
NCA, UNCA, and SUNCA. The Kondo temperature is de- 
termined as the HWHM of the Kondo peak (see inset). It is 
seen that in NCA the Kondo peak width comes out orders of 
magnitude too low. 



K 



-tt/1 



where 



r 



r 



E d + U 



(13) 



(14) 



we determine Tk as the half width of the Kondo reso- 
nance at half maximum (HWHM). 
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FIG. 6. Kondo Temperature for various parameters Ed, 
U and fixed F. Solid lines represent the exact results, Eqs. 
(13), (14). Data points are the SUNCA results determined 
from the width of the Kondo peak in the d-electron spectral 
function. 

In Fig. H the results for Tk/T obtained in this way 
for a fixed value of T = 0.05 (in units of D) and several 
values of Ed/T, as a function of I(U/T) (data points) 
are compared with the exact values Eqs. (|l3|), ( plf ) (solid 
lines). The agreement is excellent, demonstrating that 
the SUNCA provides the correct scale Tk for a wide 
range of parameters Ed and U. 
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VI. CONCLUSION 

In this paper we have proposed a conserving scheme 
to describe the Anderson impurity model at finite on- 
site repulsion U within the auxiliary particle method. 
In order to incorporate the correct value of the spin ex- 
change coupling J into the theory and, hence, to obtain 
the correct size of the low-energy scale Tk, it is neces- 
sary to treat fluctuation processes into the empty and 
into the doubly occupied intermediate state on equal 
footing at the level of bare perturbation theory. The 
simplest Luttinger-Ward functional which is completely 
symmetric in this respect consists of an infinite series of 
skeleton diagrams, corresponding to ladder-type vertex 
renormalizations in the self-energies. Although consider- 
ably more involved than the regular NCA, this approxi- 
mation, termed "symmetrized finite-C/ NCA" (SUNCA), 
is numerically tractable on a typical workstation. We 



find that SUNCA recovers the correct Kondo tempera- 
ture over a wide range of the parameters of the Ander- 
son model Ed, U , and T, while simplified approximations 
(NCA, UNCA) produce a low-energy scale typically or- 
ders of magnitude smaller than the exact value. This re- 
sult is especially relevant for a correct description of the 
low temperature properties of strongly correlated lattice 
models by a diagrammatic many-body technique, since 
in the limit of infinite dimensions these models reduce to 
a self consistent, finite-C/ single-impurity problem. Appli- 
cations of the present theory to such models are currently 
in progress. 

We are grateful to T. A. Costi and Th. Pruschke for 
useful discussions. This work was supported by DFG 
through Sonderforschungsbereich 195 (S.K., J.K., P.W.) 
and by the ESF program "Fermi liquid instabilities in 
correlated metals" (FERLIN). 



APPENDIX A: SUNCA EQUATIONS 

In this appendix we explicitly give the self-consistent 
SUNCA equations which, together with the definitions 
Eqs. (||) of the Green's functions, determine the auxil- 
iary particle self-energies. We also give the expression 
for the physical d-electron spectral function in terms of 



T a(T (uj,ft) 


= T i 


TT 




+t j 


TT 




= t j 


TT 




+t j 


TT 



where /(e) is the Fermi function, A° ca (e) — 
^ImG" (T (e)/A/'(0) the bare conduction electron density 
of states per spin, normalized to the density of states at 



the auxiliary particle propagators. 

We first define the ladder vertex functions T a , T b with 
heavy boson a and light boson b rungs, respectively, as 
shown diagrammatically in Fig. [|. These vertex func- 
tions, projected onto the physical subspace Q — 1 and 
analytically continued to real frequencies, obey the fol- 
lowing Bethe-Salpeter equations, 



(Al) 



(A2) 

I 

the Fermi level and, for concreteness, all propagators are 
to be understood as the advanced ones. The auxiliary 
particle self-energies (Fig. ^ are then given by, 



- Q)A°_ CT (e - Q)Gf- a (e)G a (e + lu - Q) 

- n)A° c _Je - ft)G/_ CT (e)G (e + lu - ft)T a _ CT ( e , fi) 

- n)A ca (fl - e)G f - a (e)G b (e + uj - fi) 

- n)A° ca (n - e)G f - a (e)G b (e + u>- 0)T b _ CT (e, O) , 



E/«rH=r J ^f(e-u)A ca (cj-e)G b (e)[l+T aa (u;,e)} 2 +T J ^f(e-co)A c _ a (e-u;)G a (e)[l + T ba (co,e)} 2 

2r 2 J ^f(e - u;)A°> - e)G b (e) J ^-f(e< - e)A°_ CT (e' - e)G f ^(e')G a (e> +cj-e) (A3) 

s 6 M = r£ / J/( e - w )^( e - w )G /ff ( e )[i + r oa ( c , w )] 

a 

J ^-f(e> - Lo)A° c _ a (e' - u)G f ^(e')G a (e' + e - w) {[1 + T ba {e, e' + e- lu)] [1 + T b _ ff (e', e' + e - w)] - 1} 



G 



£ a H=r£ f ^f( e - LJ )A° c _ (T (u J -e)G fa (e)[l + T bcr (e,u;)] 

+ T 2 E / " <")4!-> ~ *)G f *(e) x (A5) 

J ^-f(e< - a,)A°> - e')G f _ a (e')G b (e' + e - w) {[1 + r a(T (e, e' + e - «)] [1 + T a _ CT (e', e' + e - «)] - 1} . 

In order to calculate the physical impurity electron spectral function A^o- from the selfconsistently determined G a , 
Gb, Gf, it is convenient to define modified vertex functions as 





l- 


hr j 


TT 




ft)Re{G /(T (e)[l + 


T aff (e,r!)]}G a (e + W ) 


(A6) 


sU">n) = 


i - 


hr j 






0) Im{G /ff (e)[l + 


T aCT ( e ,0)]}G a (e + ^) 


(A7) 




i- 


hr j 






-e)Re{G /CT (e)[l 


+ T fcCT (e,0)]} G b (e-uj) 


(A8) 


sU",n) = 


i- 




7T 




-e) Im{G /ff (e)[l 


+ T b(7 (e,n)]}G b (e-Lo) . 


(A9) 



The impurity spectral function then reads 
1 f rffi p-@ n 

A da (u) = --Im y — _ - G /(T (fi + w) {MG b (n)]{S*_ a (u, n) 2 - n) 2 } + 2Re[G 6 (0)]S* >, fi)S£_>, 0)} 

— Im / ^ G/.^n - w) {lm[G a (tt)][S£ («, fi) 2 - ft) 2 ] + 2Re[G (fi)]S£(w, fi)5£>> J2)} 



+2 



T f dQ e~P n f de 

TT 



TT J TT f(u) 



Note that the exponentia l dive rgencies of the statistical 
factors appearing in Eq. (AlO) are compensated by the 
threshold behavior of the corresponding auxiliary parti- 
cle spectral functions ^4 M (w) = -ImG^w), fi — a,b,f 
in the integrands. For the numerical treatment, these 
divergencies can be explicitly absorbed by formulating 
the self-consistency equations (Al)-(AIO) in terms of the 



f(e - tyA c _ a (e ~ Q) lm[G b (n)G f - a {e)] Im[G /(T (0 + u)G a (e + u)] 

functions A^(uj) which are defined via 



(AlO) 



(All) 



and, hence, have no exponential divergence. We thus 
have, e.g., exp(— ^u>)A li (w) — f (ui) JLJcu) . Details of this 
representation are described in Ref.til. 
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